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RESONANCE EXPANSIONS FOR TENSOR-VALUED WAVES ON 
ASYMPTOTICALLY KERR-DE SITTER SPACES 


PETER HINTZ 


Abstract. In recent joint work with Andras Vasy [19], we analyze the low 
energy behavior of differential form-valued waves on black hole spacetimes. In 
order to deduce asymptotics and decay from this, one in addition needs high 
energy estimates for the wave operator acting on sections of the form bundle. 
The present paper provides these on perturbations of Schwarzschild-de Sitter 
spaces in all spacetime dimensions n > 4. In fact, we prove exponential decay, 
up to a finite-dimensional space of resonances, of waves valued in any finite 
rank subbundle of the tensor bundle, which in particular includes differential 
forms and symmetric tensors. As the main technical tool for working on vector 
bundles that do not have a natural positive definite inner product, we introduce 
pseudodifferential inner products , which are inner products depending on the 
position in phase space. 


1. Introduction 

We continue the analysis of (linear) aspects of the black hole stability problem 
in the spirit of earlier works by Dyatlov, Vasy and the author [10, 17, 18, 27] by 
studying linear tensor-valued, wave equations on perturbations of Schwarzschild-de 
Sitter spaces with spacetime dimension n > 4; in particular, this includes wave 
equations for differential forms and symmetric 2-tensors. In our main result, we 
establish exponential decay up to a finite-dimensional space of resonances: 

Theorem 1. Let (AI,g) denote a Kerr-de Sitter spacetime in n > 4 spacetime 
dimensions, with small angular momentum. Let £ C 71- be a subbundle of the 
bundle Tk of (covariant) rank k tensors on M, so that the tensor wave operator 
□ s = — tr V 2 acts on sections of £; for instance, one can take £ to be equal to Tk, 
symmetric rank k-tensors or differential forms of degree k. Let £1 denote a small 
neighborhood of the domain of outer communications, bounded beyond but close to 
the cosmological and the black hole horizons by spacelike boundaries, and let i* be 
a smooth time coordinate on £1. See Figure 1 for the setup. 

Then for any f £ C£°( 0,,£), the wave equation D g u = f has a unique global 
forward solution (supported in the causal future of swppf) u £ C°°(£l,£), and u 
has an asymptotic expansion 

N rrij — 1 dj 

" = EEE e - lUa n™u jm iaj m t{x) + u', 

j—1 m—0 t—1 
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where Uj m (, £ C, the resonant states aj m e, only depending on \3 g , are smooth func¬ 
tions of the spatial coordinates and <jj £ C are resonances with Imer^ > —6 (whose 
multiplicity is rrij > 1 and for which the space of resonant states has dimension dj), 
while u' £ e^ St,, L a °(n,£) is exponentially decaying, for 5 > 0 small; we measure the 
size of sections of £ by means of a t *-independent positive definite inner product. 

The same result holds true if we add any stationary 0-th order term to □, and 
one can also add stationary first order terms which are either small or subject to 
natural, but somewhat technical condition, which we explain in Remark 4.9. In fact, 
we can even work on spacetimes which merely approach a stationary perturbation 
of Schwarzschild-de Sitter space exponentially fast. See Section 2 for the form of 
the Schwarzschild-de Sitter metric and the precise assumptions on regularity and 
asymptotics of perturbations, for details on the setup, and Theorem 2.1 for the full 
statement of Theorem 1. 



Figure 1. Setup for Theorem 1 and Theorem 2.1 below. Shown 
are the black hole horizon r H + and cosmological horizon T-L , be¬ 
yond which we put an artificial spacelike hypersurface H 2 with two 
connected components. The hypersurface H 1 plays the role of a 
Cauchy hypersurface, and the forcing as well as the solution to the 
wave equation are supported in its causal future. The domain Q is 
bounded by the hypersurfaces H 1 and Hi- The ‘point at future in¬ 
finity’ in the usual Penrose diagrammatic representation is shown 
blown-up here, since the wave operator is well-behaved (namely, 
a b-operator in the sense of Melrose [23]) on the blown-up space, 
and the asymptotic information is encoded on the front face ff of 
the blow-up. 

The resonances and resonant states depend strongly on the precise form of the 
operator and which bundle one is working on. In the case of the trivial bundle, 
thus considering scalar waves, they were computed in the Kerr-de Sitter setting 
by Dyatlov [10], following work by Sa Barreto and Zworski [3] as well as Bony and 
Hafner [5]. In recent work with Vasy [19], we compute the resonances for the Hodge 
d’Alembertian on differential forms, which equals the tensor wave operator plus a 
zeroth order curvature term: We show that there is only one resonance oq = 0 in 
Imcr > 0, with multiplicity ni\ = 1, and we canonically identify the 0-resonant 
states with cohomological information of the underlying spacetime. Note however 
that [19] deals with a very general class of warped product type spacetimes with 
asymptotically hyperbolic ends, while the present paper is only concerned with 
(perturbations of) Schwarzschild-de Sitter spacetimes. We remark that in general 
one expects that D g = — tr V 2 on a bundle £ as in Theorem 1 has resonances in 
Imcr > 0, thus causing linear waves to grow exponentially in time. 
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We point out that if there are no resonances for D g (plus lower order terms) 
in Imcr > 0, thus solutions decay exponentially, we can combine Theorem 1 with 
the framework for quasilinear wave-type equations developed by the author [16] 
and in collaboration with Vasy [18] and immediately obtain the global solvability 
of quasilinear equations. This also works if there is merely a simple resonance at 
cr = 0 which is annihilated by the nonlinearity. 

The point of view from which we approach the proof of Theorem 1 was originally 
developed by Vasy [27] and extended by Vasy and the author [17]. In the context 
of scalar waves, more general and precise versions of Theorem 1 are known, see the 
references below. Thus, the main novelty is that we give a conceptually transparent 
framework that allows us to deal with tensor-valued waves on black hole spacetimes, 
where the natural inner product on the tensor bundle induced by the spacetime 
metric is not positive definite. The central motivation for the study of such waves is 
the black hole stability problem , see the lecture notes by Dafermos and Rodnianski 
[8] for details. Notice that in order to obtain energy estimates for waves, one 
needs to work with positive inner products on the tensor bundle, relative to which 
however □ is in general not well-behaved: Most severely, it is in general far from 
being symmetric at the trapped set, which prevents the use of estimates at normally 
hyperbolic trapping. In the context of black hole spacetimes, such estimates were 
pioneered by Wunsch and Zworski [29] and Dyatlov [13, 14]. On a pragmatic 
level, we show that one can conjugate □ by a suitable 0-th order pseudodifferential 
operator so as to make the conjugated operator (almost) symmetric at the trapped 
set with respect to a positive definite inner product, and one can then directly apply 
Dyatlov’s results [14]. 1 The conceptually correct point of view to accomplish this 
conjugation is that of pseudodifferential inner products , which we introduce in this 
paper. 

Roughly speaking, pseudodifferential inner products replace ordinary inner prod¬ 
ucts f{B q(u), v) \dg\ 7 where Bq is an inner product on the fibers of £, mapping £ 
into its anti-dual £ , by ‘inner products’ of the form f(B(x, D)u, v) \dg\, where 

B e is a zeroth order pseudodifferential operator mapping sections of £ into 
—% ... • • • • 
sections of £ . Thus, we gain a significant amount of flexibility, since we can allow 

the inner product to depend on the position in phase space , rather than merely on 

the position in the base: Indeed, the principal symbol b = <to(P) is an inner product 

on the vector bundle tt*£ over T*M\ 0, where it : T*M \ 0 — > M is the projection. 

One can define adjoints of operators P £ 4/ m (M ,£) (e.g. P = D g ), acting 
on sections of £, relative to a pseudodifferential inner product B , denoted P* B , 
which are well-defined modulo smoothing operators. Moreover, there is an invariant 
symbolic calculus involving the subprincipal operator S su b(P) 7 which is a first order 
differential operator on T*M\0 acting on sections of n*£ that invariantly encodes 
the subprincipal part of P, for computing principal symbols of commutators and 
imaginary parts of such operators. In the case that P is principally scalar and real, 
the principal symbol of P — P* B £ 4' m_1 (M, £) then vanishes in some conic subset 
of phase space T*M \ 0 if and only if S' su b(P) — Ssub{P)* b does, which in turn 
can be reinterpreted as saying that the principal symbol of QPQ~ X — (QPQ~ 1 )* B ° 
vanishes there, where Bq is an ordinary inner product on £, and Q £ 4 >°{M,£) is 
a suitably chosen elliptic operator. In the case considered in Theorem 1 then, it 

1 1n other words, we reduce the high frequency analysis of tensor-valued waves to an essentially 
scalar problem. 
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turns out that the subprincipal operator of D g on tensors, decomposed into parts 
acting on tangential and normal tensors according to the product decompositions 
M = R t x X x and X = (r_,r+) x S” -2 , at the trapped set equals the derivative 
along the Hamilton vector field Hq, G the dual metric function, plus a nilpotent 
zeroth order term. This then enables one to choose a positive definite inner product 
b on 7 t*£ relative to which 5 , su b(Dg) is arbitrarily close to being symmetric at the 
trapped set; thus with B = b(x,D), the operator D g is arbitrarily close to being 
symmetric with respect to the pseudodifferential inner product B. Hence, one can 
indeed appeal to Dyatlov’s results on spectral gaps by considering a conjugate of 
□ g , which is the central ingredient in the proof of Theorem 1. 

1.1. Related work. The study of non-scalar waves on black hole backgrounds has 
focused primarily on Maxwell’s equations: Sterbenz and Tataru [25] showed local 
energy decay for Maxwell’s equations on a class of spherically symmetric asymptoti¬ 
cally flat spacetimes including Schwarzschild. Blue [4] established conformal energy 
and pointwise decay estimates in the exterior of the Schwarzschild black hole; An- 
dersson and Blue [1] proved similar estimates on slowly rotating Kerr spacetimes. 
These followed earlier results for Schwarzschild by Inglese and Nicolo [22] on energy 
and pointwise bounds for integer spin fields in the far exterior of the Schwarzschild 
black hole, and by Bachelot [2], who proved scattering for electromagnetic pertur¬ 
bations. Finster, Kamran, Smoller and Yau [15] proved local pointwise decay for 
Dirac waves on Kerr. There are further works which in particular establish bounds 
for certain components of the Maxwell field, see Donniger, Schlag and Softer [9] and 
Whiting [28]. Dafermos [6], [7] studied the non-linear Einstein-Maxwell-scalar field 
system under the assumption of spherical symmetry. 

The framework in which we describe resonances was introduced by Vasy [27]. 
In the scalar setting, this can directly be combined with estimates at normally 
hyperbolic trapping [13, 14, 29] to obtain resonance expansions for scalar waves. 
On exact Kerr-de Sitter space, Dyatlov proved a significant strengthening of this in 
[11], obtaining a full resonance expansion for scalar waves, improving on the result 
of Bony and Hafner [5] in the Schwarzschild-de Sitter setting, which in turn followed 
Sa Barreto and Zworski [3] . Vasy [26] proved the meromorphic continuation of the 
resolvent of the Laplacian on differential forms on asymptotically hyperbolic spaces, 
and the fact that the underlying analysis of [27] works on sections of vector bundles 
just as it does on functions is fundamental for the present paper. 

1.2. Structure of the paper. In Section 2, we recall the Schwarzschild-de Sitter 
metric and its extension past the horizons, put it into the framework of [17, 27] for 
the study of asymptotics of waves, and establish the normally hyperbolic nature of 
its trapping. We proceed to sketch the proof of Theorem 1, leaving the discussion of 
high energy estimates at the trapped set to the subsequent sections, which comprise 
the central part of the paper: We introduce pseudodifferential inner products on 
vector bundles in full generality in Section 3, and we use the theory developed 
there in Section 4 to study pseudodifferential inner products for wave operators on 
tensor bundles, uncovering the nilpotent nature of the subprincipal operator of □ 
on Schwarzschild-de Sitter space at the trapping in Section 4.2 and thereby finishing 
the proof of Theorem 1. 

Acknowledgments. I am very grateful to Semyon Dyatlov, Andras Vasy and 
Alexandr Zamorzaev for many useful discussions. 
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2. Detailed setup and proof of the main theorem 


We recall the form of the ?x-dimensional Schwarzschild-de Sitter metric, n > 4: 
We equip M = I t x X, X = (r_,r+) r x S™ -2 , with r± dehned below, with the 
metric 

go = pdt 2 — (/x _1 dr 2 + r 2 dco 2 ), ( 2 . 1 ) 

where dui 2 is the round metric on the sphere § n_2 , and // = 1 — — Ar 2 , 

A = („_ 2 )pi.-i) ’ w ith M, > 0 the black hole mass and A > 0 the cosmological 
constant A. The assumption 


Ml A n “ 3 < ^ 

(n 


3)n — 3 

l ) 71 " 1 


( 2 . 2 ) 


guarantees that /i has two unique positive roots 0 < r_ < r+. Indeed, let pi = 
r~ 2 p = r~ 2 - 2M.r 1 ~ n - A. Then ft' = - 2r~ n (r n " 3 - (n - 1 )M.) has a unique 
positive root r p — [(n — l)-/!/.] 1 ^" -3 ), j2'(r) > 0 for r £ (0,r p ) and ju'(r) < 0 
for r > r p ; moreover, Jl(r) < 0 for r > 0 small and Jl(r) -A —A < 0 as r —► oo, 
thus the existence of the roots 0 < r_ < r + of Jl is equivalent to the requirement 
Jl{ r p ) = r~ 2 — A > 0, which is equivalent to (2.2). 

Define a = /W 2 , thus da = \[>!a~ l dr, and let 

^ ±(r):=T /T(r) (2 ' 3) 


near r±, so /3±(r±) > 0 there. Then the metric go can be written as 

go = a 2 dt 2 — h, h = a~ 2 dr 2 + r 2 duj 2 = da 2 + r 2 du 2 , 

We introduce a new time variable t* = t — F(a), with d a F = — a - 1 /3± near r = r±. 
Then 

go = At dt 2 — (3± dt* d/x — r 2 dui 2 

near r = r±, which extends as a non-degenerate Lorentzian metric to a neighbor¬ 
hood M = K ti> x X of M, where X = (r_ — 25, r+ + 2d) x S n-2 . We will consider 
the Cauchy problem for the tensor wave equation in the domain f 1 C M, 

= [0, oo)t„ x [r_ — 5, r + + 5} r x S n ~ 2 . 


Thus, H bounded by the Cauchy surface Fd i = {t* = 0}, which is spacelike, and 
by the hypersurface H 2 = U±{ r = r±± <5}, which has two spacelike components, 
one lying beyond the black hole (r_) and the other beyond the cosmological (r+) 
horizon; see Figure 1. 

For the purpose of analysis on spacetimes close to (but not necessarily asymptot¬ 
ically equal to!) Schwarzschild-de Sitter space, we encode the uniform (asymptot¬ 
ically stationary) structure of the spacetime by working on a compactified model, 
which puts the problem into the setting of Melrose’s b-analysis, see [23]: Define 
r := e _t *, and bordify M to a manifold M with boundary by adding r = 0 as 
the boundary at future infinity and declaring r to be a smooth boundary defin¬ 
ing function. The metric go becomes a smooth Lorentzian b-metric on M: If 

— . . . . . 2 
dxi denotes coordinate differentials on X, then g 0 is a linear combination of a \ , 

— <g) dx^ + dxi ( 8 > — and dxi 8 > dxj with coefficients which are smooth on M, and 
go, written in such coordinates, is a non-degenerate matrix (with Lorentzian sig¬ 
nature) up to and including r = 0. Invariantly, we have the Lie algebra Vb(M) 
of b-vector fields, which are the vector fields tangent to the boundary, spanned by 
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rd T = —dt, and d Xi ; elements of Vb(M) are sections of a natural vector bundle 
b TM, the b-tangent bundle, and we have the dual bundle h T*M , spanned by ^ 
and dXi. Thus, g is a smooth non-degenerate section of the symmetric second tensor 
power S 2h T*Wl. 

Now, given a complex vector bundle £ —> M of finite rank, equip it with an arbi¬ 
trary Hermitian inner product and any smooth b-connection, which gives a notion 
of differentiating sections of £ along b-vector fields; over £1 (which has compact 
closure in M), all choices of inner products are equivalent. We can then define the 
b-Sobolev space H^(£l,£) for s £ Z>o to consist of all sections of £ over Id which 
are square integrable (with respect to the volume density \dg\ induced by the metric 
g ) together with all of its b-derivatives up to order s, and extend this to all s S M 
by duality and interpolation, or via the use of b-pseudodifferential operators. For 
the forward problem for the wave equation, we work on spaces of functions which 
vanish in the past of Hi and which extend across H 2 - Thus, we work with the space 
H^(£l,£)*~ of distributions u £ H^(£l,£) which are extendible distributions at H 2 
and supported distributions at H i, i.e. they are restrictions to f 2 of distributions on 
M which are supported in t* > 0. See Hormander [21, Appendix B] for details. We 
also have weighted b-Sobolev spaces H^’ r (£l,£) = r r H^(£l,£), likewise for spaces 
of supported/extendible distributions. Note that the b-Sobolev spaces H ^ are inde¬ 
pendent of the choice of boundary defining function r in that the choice t' = ar 7 , 
a = a(x) smooth, 7 > 0, while it changes the smooth structure of M, yields the 
same spaces with equivalent norms. The asymptotic behavior of waves will be 
encoded on the boundary doo£1 at future infinity of £1, that is, on 

doc£1 = {t = 0} x [r_ — 6, r + + <5] r x § n ~ 2 , 

which is a smooth manifold with boundary. Similarly to the above definitions, 
we can define Sobolev spaces (including semiclassical versions of these) with sup¬ 
ported/extendible character at the boundary. 

Suppose g is a Lorentzian b-metric such that for some smooth Lorentzian b- 
metric g', we have g — g' £ H £°'’ (Q, S' 2h T*M) for some r > 0 . 2 Changing g' so as 
to make it invariant under time translations does not affect this condition, so let us 
assume g' is i*-invariant. We will consider the wave operator D g acting on sections 
of the bundle Tk of covariant tensors of rank k over O. We assume that g' and 
go are close (in the C k sense for sufficiently high k), 3 so that the dynamical and 
geometric structure of g is close to that of go' 4 Most importantly, the nature of the 
trapping for g' (and thus for g) is still normally hyperbolic , 5 and the subprincipal 
operator (see Section 3.3) of □, at the trapped set, while not necessarily having the 
nilpotent structure alluded to in the introduction and explained in Section 4.2, has 


2 By the discussion of b-Sobolev spaces above, this condition on g is invariant, i.e. independent 
of the specific choice of the boundary defining function e~ t " of the spacetime at future infinity. 

5 In other words, the metric g is exponentially approaching a stationary metric close to the 
Schwarzschild-de Sitter metric, so for instance perturbations (within this setting) of Kerr-de Sitter 
spaces are allowed. 

^See [17, §3] and [18, §5] for details. 

5 We will show the r-normal hyperbolicity for every r of the trapping for Schwarzschild-de 
Sitter in all spacetime dimensions below, and r-normal hyperbolicity (for large, but finite r) is 
structurally stable under perturbations of the metric, see Dyatlov [13] and Hirsch, Shub and Pugh 
[ 20 ]. 
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small imaginary part relative to (the symbol of) a pseudodifferential inner product 
on 7fc. We then have: 

Theorem 2.1. In the above notation, if g' is sufficiently close to the Schwarzschild- 
de Sitter metric go, then there exist So € R and 6 > 0 as well as a finite set 
{<jj: j = 1,..., TV} C C, Irn <jj > —5, integers mj > 1 and dj > 1, and smooth 
functions aj m f, £ C°° (docfl), l<j<N,0<m< mj — 1, 1 < I < dj, such that the 
following holds: The equation 

a g u = f, f £H^ s s > s 0 , (2.4) 

has a unique solution u £ Hff 00 ' °°(^j7fc)*’ - > which has an asymptotic expansion 

N rrij — 1 dj 

u = X(t) ^ Y I lo g + u', 

j—1 m =0 £=1 

where \ is a cutoff function, i.e. x(r) = 1 near t = 0 andx( T ) = 0 near the Cauchy 
surface Hi, and Uj m e £ C, while the remainder term is u' £ H(j * S (VL,Tk)*'~ ■ 

The same result holds true if we restrict to a subbundle ofTk which is preserved 
by the action ofO, for instance the degree k form bundle, or the symmetric rank k 
tensor bundle. 

If V £ C°°(M, End(7I-)) + End(7I-)), r > 0, is a smooth (conormal) 

End(7fc )-valued potential (without restriction on its size), the analogous result holds 
for D g replaced by D g + V. We may even change Dg by adding a first order b- 
differential operator L acting on T( with coefficients which are elements of C°° + 
H^ ,r , provided either the coefficients of L are small, or the subprincipal operator 
ofOg + L is sufficiently close to being symmetric with respect to a pseudodifferential 
inner product on Tk, see Remark 4.9. 

The numbers <jj are called resonances or quasinormal modes, and the functions 
ajmt resonant states. They have been computed in various special cases; see the 
discussion in the introduction for references. The threshold regularity s o is related 
to the dynamics of the flow of the Hamiltonian vector field of the dual metric 
function G (i.e. G(x,f) = |£|g, ), with G the dual metric of g) near the horizons 
which are generalized radial sets , see [17, Proposition 2.1]. Thus, so can easily be 
made explicit, but this is not the point of the present paper. 

The proof of Theorem 2.1 proceeds in the same way as the proof of [17, Theo¬ 
rem 2.20] in the scalar setting, so we shall be brief: Denote by N(O g ) the normal 
operator of D g : We freeze the coefficients of Dg £ Diff^(D,7fc) at dooll and thus 
obtain a dilation-invariant operator N(D g ), with D g — N(D g ) being an operator 
whose coefficients decay exponentially (in f*) by assumption on the structure of g. 
Denote by □g(<r) £ Diff 2 (9 00 H, Tk) the Mellin transformed normal operator family, 
depending holomorphically on a £ C, which we obtain from N(D g ) by replacing 
D t * by — cr. 6 Once we show high energy estimates for □(cr)~ 1 , which are polynomial 
bounds on its operator norm between suitable Sobolev spaces' as |Recr| —> oo in 
Im a > — 5 , we can use a contour shifting argument to iteratively improve on the 

^Changing the boundary defining function r to a(.r)r 7 , we can express the normal operator 

with respect to the new defining function in terms of the normal operator with respect to r, 

namely it equals a(x)~ 1 £ig('yo')a(x). 

7 These are semiclassical Sobolev spaces with extendible character at the boundary of 
see in particular [27] and the proof of [17, Theorem 2.20]. 
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decay of u, picking up contributions of the poles of □(er)~ 1 which give rise to the 
resonance expansion. 8 The fact that the remainder term u' has the same regularity 
as the forcing term /, thus u' loses 2 derivatives relative to the elliptic gain of 
2 derivatives, comes from the high energy estimate losing a power of 2, see [18, 
Theorem 5.5], which in turn is caused by the same loss for high energy estimates 
at normally hyperbolic trapping, see [14, Theorem 1], or [18, Theorem 4.5] for a 
microlocalized version of Dyatlov’s estimate. 

Thus, the crucial point is to obtain high energy estimates at the trapped set for 
the operator □ acting on Tk in liner > —S. Dyatlov’s result [14, Theorem 1] (see 
also the discussion preceding [18, Theorem 5.5]) shows that a sufficient condition 
for these to hold is 

|0-| -1 CT b ,l (^7 (□-□*))< Pmin/2 (2.5) 

at the trapped set T, 9 where i/ m ; n is the minimal normal expansion rate of the 
Hamilton flow at the trapping, see [14] and the computation below. Here, the 
adjoint is taken with respect to a positive definite inner product on Tk', note that 
the inner product induced by g, with respect to which □ is of course symmetric, is 
not positive definite, except when k = 0, i.e. for the scalar wave equation. Since g 
is close to the Schwarzschild-de Sitter metric, it suffices (by the dynamical stability 
of the trapping) to obtain such a bound for the Schwarzschild-de Sitter metric go- 
While this bound is impossible to obtain directly for the full range of Schwarzschild- 
de Sitter spacetimes, we show in Section 4.2 how it can be obtained if we use 
pseudodifferential products. Prosaically, this means that we consider a conjugated 
operator P := QOQ~ , where Q £ \t°(M, Tk) is elliptic with parametrix Q~ , and 
for any e > 0, we can arrange |cr| _1 crb,i(7^(P — P*)) < e (with the adjoint taken 
relative to an ordinary positive definite inner product on Tk), thus (2.5) holds for 
□ replaced by P; we will prove this in Theorem 4.8. Hence [14, Theorem 1] applies 
to P, establishing a spectral gap; indeed, by the remark following [14, Theorem 1], 
Dyatlov’s result applies for operators on bundles as well, as soon as one establishes 
(2.5). Arranging (2.5) in a natural fashion lies at the heart of Sections 3 and 4. 

It remains to establish the r-normal hyperbolicity for all r for the Schwarzschild- 
de Sitter metric. The dynamics at the trapping only depend on properties of the 
(scalar!) principal symbol go of □. For easier comparison with [12, 27, 29], we 
consider the operator P = —r 2 D instead. We take the Fourier transform in —t 
and rescale to a semiclassical operator on X (this amounts to multiplying P by h 2 , 
giving a second order semiclassical differential operator P^, with h = |oj _1 , and 
we then define z = her). Introducing coordinates on T*X by writing 1-forms as 
£ dr + gcLo, and letting 

A r = r 2 n = r 2 ( 1 - Ar 2 ) - 2 M.r 5 * * ~ n , 


®As shown by Vasy [27, §7], these estimates in Imcr 0 are automatic if the boundary defining 

function of future infinity is timelike; our choice does not satisfy this, but changing t* by a smooth 

function of the spatial variables, this can easily be arranged, see [27, §6], and in fact we can arrange 

f* = t away from the black hole and cosmological horizons. 

9 We work in the b-setting, which via the Mellin transform is equivalent (on the normal operator 
level, which is all that matters) to the semiclassical setting considered in Dyatlov’s work, see the 
discussion in [18, §5]. 
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Ph has semiclassical principal symbol 

p=A r e-^-z 2 +\ v \ 2 , 

and correspondingly the Hamilton vector field is 

H p — 2A r £d r — ^ d r A r l; 2 — + H\ v \2 

We work with real z, hence z = ±1. First, we locate the trapped set: If H p r = 
2A r £ = 0, then £ = 0, in which case H 2 r = 2A r H p £ = 2A r d r (r 4 /A r )z 2 . Recall 
the definition of the function J1 = \xlr 2 = A r /r , then we can rewrite this as 
H 2 r = —2A r Jl~ 2 (d r Jl)z 2 . We have already seen that d r Jl has a single root r p £ 
(r_,r_(_), and (r — r p )d r jl < 0 for r ^ r p . Therefore, H p r = 0 implies (still assuming 
H p r = 0) r = r p . We rephrase this to show that the only trapping occurs in the 
cotangent bundle over r = r p : Let F(r) = (r — r p ) 2 , then H p F = 2 (r — r p )F[ p r 
and H p F = 2(H p r) 2 + 2(r — r p )H 2 r. Thus, if H p F = 0, then either r = r p , 
in which case H 2 F = 2(H p r) 2 > 0 unless H p r = 0, or H p r = 0, in which case 
H p F = 2{r — r p )H 2 r > 0 unless r = r p . So H p F = 0,p = 0 implies either H p F > 0 
or r = r p , H p r = 0, i.e. 

4 

(r,w;£,7?) <E Tft := |(r p , w; 0, ?y): ^-z 2 = |r?| 2 |, 


so T h is the only trapping in T*X, and F is an escape function. We compute the 
linearization of the H p -flow at T^ in the normal coordinates r — r p and £, to wit 


r ( r ~ r p\ _ ( 0 

P V ^ ~ \2(n-3)r- 4 (fi\ r=rp )- 2 z 2 

+ 0(\r-rp\ 2 + \tf), 



where we used d rr lAr=r p = —2 (n — 3)r p 4 , which gives <9 r /i = —2 (n — 3)r p 4 (r — 
r p ) + 0(\r — r p | 2 ). The eigenvalues of the linearization are therefore 


±2 TV 


n — 1 


1/2 


i n-l „ 2 \ 

4 T P A , 


which reduces to the expression given in [27, p. 85] in the case n = 4, where 
r p = 3 M, = | r s with r s = 2M t , and A = A/3. In particular, the minimal 
expansion rate for the semiclassical rescaling of □ at the trapping T^ is 

/ i \ 1/2 

- 11 ' > 0. 


A ni t i — 2r/, 


1 n ~ 1 -2 \ 
A F=3 r P A , 


The expansion rate of the flow within the trapped set is 0 by spherical symme¬ 
try; note that integral curves of H p on T[ are simply unit speed geodesics of the 
round unit sphere § n ~ 2 . This shows the normal hyperbolicity (in fact, r-normal 
hyperbolicity for every r) of the trapping and finishes the proof of Theorem 2.1. 

For later reference, we note that the spacetime trapped set, i.e. the set of points 
in phase space that never escape through either horizon along the Hamilton flow, 
is given by 

T = {( t,r = r p , cv; cr, £ = 0 ,r?): ct 2 = ^M 2 }, 

L ,vk'(r p )=0. 


( 2 . 6 ) 
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3. Pseudodifferential inner products 

We now develop a general theory of pseudodifferential inner products, which we 
apply to the setting of Theorem 2.1 in Section 4. 

We work on a complex rank N vector bundle £ over the smooth compact 71- 
dimensional manifold A' without boundary. We will define pseudodifferential inner 
products on £, which are inner products depending on the position in phase space 
T* A, rather than merely the position in the base X. As indicated in the intro¬ 
duction, we achieve this by replacing ordinary inner products by pseudodifferential 
operators whose symbols are inner products on the bundle n*£ —> T*X \ 0, where 
7 r: T*X \ 0 —> X is the projection. 

3.1. Notation. Let V be a complex TV-dimensional vector space. We denote by V 
the complex conjugate of V, i.e. V = V as sets, and the identity map t: V —t V is 
antilinear, so l(Xv) = A l(v) for v £ V, A £ C, which defines the linear structure 
on V . 10 A Hermitian inner product H on V is thus a linear map H: V (g> V —> C 
such that H(u,l(v)) = H(v,l(u)) for u,v £ V, and H(u,l(u)) > 0 for all non- 
zero u € V. This can be rephrased this in terms of the linear map B: V —> V 
defined by B{u) = H(u , •) and the natural dual pairing of V with V, namely 
( Bu,l(v )) = ( Bv,l(u )), and ( Bu,l(u )) > 0 for u £ V non-zero. 

A map A: V —> V has a transpose A T : V —> V *, which satisfies ( Au,l(v )) = 
(u,A t l(v)) for all u,v £ V, and an adjoint A*: V —> V satisfying ( Au,l(v )) = 
(A*v, l(u)). Concretely, defining the antilinear map 

j: V* ->• V*, (, j(£),i(v )) = (£,v), 

we have A* = jA T L. The symmetry of a Hermitian inner product B as above 
is simply expressed by B = B*. Similarly, a map P: V —> V has a transpose 
P T : V* —> V* and an adjoint P*: V* —> V* defined by (I,i(Pv)) = ( P*£,t(v )) 
for J £ V and v £ V, and one easily finds P* = jP T j~ 1 . We point out that 
the definitions of adjoints of maps A: V —> V and P: V —> V are compatible in 
the sense that ( AP )* = P*A*. Furthermore, if B: V —l V is a Hermitian inner 
product and Q: V —> V is invertible, then B\ = Q*BQ defines another Hermitian 
inner product, (Biu,l(v)) = (BQu, i(Qv)). 

Now, given an inner product B on V and any map P: V —> V, the adjoint P* B 
of P with respect to B is the unique map P* B : V —> V such that ( BPu , t( v)) = 
(B 71 , l(P* b v)) for all it, v £ V. We find a formula for P* B by computing 

(BPu, i(v)) = (B*(B*)- 1 P*B*v,l{u)) = (Bu , l((B PB- 1 )*v)) , 

i.e. P* B = ( BPB -1 )* = B~ 1 P*B. The self-adjointness of P with respect to B is 
thus expressed by the equality P = B~ 1 P*B. 

If £ is a complex rank N vector bundle, we can similarly define the complex 

conjugate bundle £ as well as adjoints of vector bundle maps £ —> £ and £ —> 
—♦ 

£ . We can also define adjoints of pseudodifferential operators mapping between 
these bundles: For convenience, we remove the dependence of adjoints on a volume 
density on X by tensoring all bundles with the half-density bundle fls over A', and 


1 (l We prefer to write /.(?;) rather than v to prevent possible confusion with taking complex 
conjugates in complexifications of real vector spaces. 
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we have a natural pairing 

(£* <8> x (£ 0 fl® ) x 9 (?, l{v)) ha (I, i(v)) £ f2*, x £ X, 

likewise for the complex conjugate of £. Thus, an operator A G \k m (X, £® fia 5 £ 0 
ffs) has an adjoint A* £ ^f m (X,£ 0 Qs, £* <g> ) defined by 

[ (A*u,l(v)}= ( { Av,l{u )), 

Jx Jx 

with principal symbol cr m (A*) = cr m (A)* £ S m (T*X \ 0, n* Hom(f, £*)), and like¬ 
wise P £ i S> rn (X,£ 0 fi3) has an adjoint P* £ H/ m (X,£* 0 ids) with a m (P*) = 
<T m (P)*. 


3.2. Definition of pseudodifferential inner products; adjoints. We work 
with classical, i.e. one-step polyhomogeneous, symbols and operators, and denote 
by SP (T*X \ 0) symbols which are homogeneous of degree m with respect to 
dilations in the fibers of T*X \ 0. 

Definition 3.1. A pseudodifferential inner product (or tH-inner product) on the 
vector bundle £ — > X is a pseudodifferential operator B £ \I/ 0 (A';£ 0 ; £ 0 ffs) 

satisfying B = B* , and such that moreover the principal symbol cr°(B) = b £ 
^hom(^bY \ 0; 7T* Hom(£,£*)) of B satisfies 

(b(x,£)u,L(u)) > 0 (3.1) 

for all non-zero u £ £ x , where x £ X, £ £ T*X \ 0. If the context is clear, we will 
also call the sesquilinear pairing 

C°°(X, £ 0 fD) x C°°(X, £ 0 fl 5 ) 9 (u,v) ha f ( B(x,D)u,t(v )) 

Jx 

the pseudodifferential inner product associated with B. 

In particular, the principal symbol b of B is a Hermitian inner product on tt*£. 
Conversely, for any b £ S® om (T*X \ 0; n* Hom(£,£*)) satisfying b = b* and (3.1), 
there exists a 'h-inner product B with <r°(B ) = 6; indeed, simply take B to be any 
quantization of b and put B = \{B + B*). 

Remark 3.2. While we will develop the theory of \I/-inner products only in the 
standard calculus on a closed manifold, everything works mutatis mutandis in other 
settings as well. Thus, in the b-calculus of Melrose [23], Tb-inner products on a 
manifold with boundary are defined similarly to ik-inner products, except that 
adjoints are defined on the space C°° of functions vanishing to infinite order at the 
boundary, and the space of ‘trivial’, smoothing operators is now dt^ 00 , likewise for 
the scattering calculus [24], replacing ‘b’ by ‘sc’. In the semiclassical calculus on a 
closed manifold, adjoints are again defined on C °°, but the space of ‘trivial’ operators 
is now h 00 ^^ 00 , and suitable factors of h need to be put in for computations 
involving subprincipal symbols. 

We next discuss adjoints of dlDOs relative to ik-inner products. 

Definition 3.3. Let B be a T-inner product, and let P £ \k m (A', £ 0 iAi) , then 
P* B £ \l/ m (X, £ 0 f2s) is called an adjoint of P with respect to B if there exists an 
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operator R G ^ 00 (X, £ (g) 2 ? £* <g> 2) such that 

J (BPu,l(v)) = J ( Bu,l{P* b v )) + J (Ru, l(v)) (3.2) 

for all u,v GC°°(X,£0f25). 

Remark 3.4. This definition and the following lemma have straightforward general¬ 
izations to the case that P maps section of £ into sections of another vector bundle 
T, provided a (\I/-)inner product on T is given. 

Lemma 3.5. In the notation of Definition 3.3, the adjoint of P with respect to B 
exists and is uniquely determined modulo T” 00 (X, £). In fact, P = ( BPB ~)*, 
where B~ is a parametrix for B. Moreover, ( p* B y B = p modulo £0f22). 

In particular, Im B P = ^(P — P* B ) is self-adjoint with respect to B (i.e. its own 
adjoint modulo T -00 ). 

Proof. Let B~ be a parametrix of B and put Rl = I — B~B G ’F~ 00 (X ,£ 0 12 i). 
Then 

J (BPu,t(v)) = J ( BPB~Bu,l(v )) + ( BPR l u,l{v )), 

hence (3.2) holds with P* B = ( BPB ~)* and i? = BPRl- To show the uniqueness 
of P* B modulo smoothing operators, suppose that P is another adjoint of P with 
respect to B, with error term R (i.e. (3.2) holds with P* B and R replaced by P 
and R). Then 

J (B(P* B -P)v,l(u)) = j(Bu , l((P* b - P)v)) = J{(R — R)u , t («)> 

= J {{R- R)*v,l(u)) 

for u,v G C°°(X,£®ni), so B(P* b -P) = (R-R)* G T-°°(X, £ (g) 12 $,£* <8>f2s), 
and the ellipticity of B implies P* B — P G \f ,-00 (A’’ ) £ (g> 12$), as claimed. 

Since B is self-adjoint, we can assume that B~ is self-adjoint by replacing it by 
\(B~ + (B~)*) (which changes B~ by an operator in T^ 00 ). Then the second 
claim follows from 

(p*By*B = ^bP* b B~)* = B~ BPB~ B = P 
modulo £(8)125). □ 

Lemma 3.6. Suppose P G 4/ m (X, £ ® Sli) is self-adjoint with respect to B. Then 
its principal symbol p is self-adjoint with respect to b = a 0 (B), i.e. 

(b(x,£)p(x,£)u,i(v)) = (b(x,^)u,L(p(x,^)v)), x G X,£ G T x X,u,v G £ x . 

Proof. The hypothesis on P means ( BPB ~)* = P modulo 4/ _00 ; thus on the level 
of principal symbols, p = b~ 1 p*b = p * b , which proves the claim. □ 

We now specialize to the case that P G 4 ,m (X, £®f2s) has a real, scalar principal 
symbol. Fix a coordinate system of X and a local trivialization of £, then the full 
symbol of P is a sum of homogeneous symbols p ~ p m + p m -i + • •with pj 
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homogeneous of degree j and valued in complex N x N matrices. Recall from [21, 
§18] that the subprincipal symbol 

a svh {P)=p m - i(x,Z) - e SZ^(T*X\ 0,C JVxJV ) (3.3) 

j 

is well-dehned under changes of coordinates; however, it does depend on the choice 
of local trivialization of £. We compute the principal symbol of 

Im B P ■- — (P~P* B ) 

2 i ' 

for such P in a local trivialization of £; we will give an invariant formulation in 
Proposition 3.10 below. 

Lemma 3.7. Let P £ g f li) be a principally real and scalar, and let 

B = b(x,D) be a T -inner product on £. Then Im B P £ '^ m ~ 1 {X,£ g fis) has the 
principal symbol 


P) = Im b a sub (P) + -b^H^b), 


(3.4) 


where Im b cr su b(P) = ^(a'sub(-P) — cr su b(P)* b ). Here, we interpret b and a S ub(P) 
as N x N matrices of scalar-valued symbols using a local frame of £ and the corre¬ 
sponding dual frame of £ , and the action of H p is component-wise. 

Proof. We compute in a local coordinate system over which £ and £ are trivi- 
• , # ^ 
alized by a choice of N linearly independent sections ei,...,ejV) and £* and £ 

are trivialized by the dual sections e\,...,e* N £ £* satisfying e* (e.j) = 6ij, ex¬ 
tended linearly as linear functionals on £, resp. on £, in the case of £*, resp. £ . 
We trivialize fusing the section \dx\z. Let = (b(x,tf)ej,i(ei)), then 

b(x, £) = ( bij(x , t;))ij=i,...,N, a linear map from the fibers of £ to the fibers of £ , is 
the symbol of B in local coordinates: If u = ]>W u 3 e.j \dx\ ? and v = JN Vjej\dx\ s, 
we have 

(b(x,£)u,i(v)) ='^2b ij (x,^)u j ■vl\dx \, 

*3 

thus 


Note that b(x,t;) is a Hermitian matrix, i.e. bij(x,£) = bji(x,tf), and in fact B = 
b(x , D ) is self-adjoint (with respect to the standard Hermitian inner product on 
C^). The adjoint of P = p(x,D), which in local coordinates is simply an TV x N 
matrix of scalar TDOs, with respect to B is the operator P = p(x, D) such that 


Her 


J b(x, D)p(x, D)u ■ v dx = J b(x, D)u ■ p(x, D)v dx + J Ru-vdx, 

Let B~ := b~(x,D) be a parametrix for b{x,D), in particular b~ (x ,£) = ^(a:,^)^ 1 
modulo S' -1 ; we may assume B~{x,D)* = B~{x,D). We then have 

p(x, D) = b~(x, D)p(x , D)*b(x, D) 
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by Lemma 3.5. Write p(x, £) = p m {x, £) +Pm-i{x, £) + ..., then the full symbol of 
P — P = B~ {BP — P*B) (where P* is the adjoint of P with respect to the standard 
Hermitian inner product on C N ) is given, modulo S m ~ 2 , by 

b~ X (bp m • E d^bd Xj p m + bp m— 1 

3 

~ P ™ b ~) Yl^ dx j^P*m) b - \ E d ^P*m d xjb- P* m --ib) 

3 3 

= i^Pm— 1 'y ' 9xj ZjPm^ — b (^Pm— 1 y ' ^Xj^jPrn'j b + ib H p m (5), 
j j 

where we used that p m is scalar and real. The claim follows. □ 

3.3. Invariant formalism for subprincipal symbols of operators acting on 
bundles. We continue to denote by P £ \E' m (A, £ (§) 12 2) a principally scalar 'I'DO 
acting on the vector bundle £, with principal symbol p. 11 We will show how to 
modify the definition (3.3) of the subprincipal symbol of P, expressed in terms 
of a local trivialization of £, in an invariant fashion, i.e. in a way that is both 
independent of the choice of local trivialization and of local coordinates on X. This 
provides a completely invariant formulation of Lemma 3.7. 

Let U C X be an open subset over which £ is trivial, and pick a frame e(x) = 
{ei(x ),..., ejv(a:)} trivializing £ over U. Let us write P e for P in the frame e, i.e. 
P e = (Pj k )j,k=i,...,N is the TV x N matrix of operators PJ k £ \E ,m (Z7, fis) defined by 

F (yZk u k{x)e k (x)) = J2jk P jk( u k) e j { x )> u k £ Then of ub (P) as defined 

in (3.3), with the superscript making the choice of frame explicit, is simply an 
N x N matrix of scalar symbols: 

cr sub(P) = ( a sub(Pjk))j,k—l,...,N- 

We will consider the effect of a change of frame on the subprincipal symbol (3.3). 
Thus, let C £ C°°(U, End(£)) be a change of frame, i.e. C(x) is invertible for all 
x £ A'. Then ej(x) = C(x)eb(x) defines another frame e'(x) = {e^(cc),... ,e' N (x)} 
of £ over U. One easily computes 

oi^C-'PC) = (Ce'r'ai^P)^’ -i(C e ')- 1 H p (C e '), 

with H p interpreted as the diagonal TV x N matrix Inxn^p of first order differential 
operators, and C e is the matrix of C in the frame e!. Now note that ( C~ 1 PC) e = 
P e and ( C e ) -1 H p (C e ) = ( C e ) _1 H p C e — H p ; thus, we obtain 

<ub (P) - = {C e ')-\cri h {P) - iH p )C e ' (3.5) 

Thus, viewing cr® ub (P) — i H p as the TV x TV matrix (in the frame e') of a differential 
operator acting on C°°(T*X \ 0,7 t*£), the right hand side of (3.5) is the matrix 
of the same differential operator, but expressed in the frame e. Notice that the 
principal symbol p of P as a scalar, i.e. diagonal, TV x TV matrix of symbols, is 
well-defined independently of the choice of frame. To summarize: 


^The discussion until Proposition 3.8 works for principally non-scalar operators as well with 
mostly notational changes. 
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Definition 3.8. For P £ *F m (X, £ (g) £7i) with scalar principal symbol p, there is a 
well-defined subprincipal operator S sn \,(P) £ Diff 1 (T*X \0,7r*f), homogeneous of 
degree m — 1 with respect to dilations in the fibers of T*X \ 0, defined as follows: 
If (ei(a;),... ,ejv(x)} is a local frame of £, define the operators Pjk £ 'F m (X, fii) 
by p (J2k u k{x)e k (x)) = J2jk P jk(uk)ej(x), u k G C°°(X,£l5). Then 

*5sub (P) ^ q k (,X, ^ ’ 1 U.sn h ( l\jk ) Uk ) Cj j ^ ) (H pQk'jCk- 

k jk k 

In shorthand notation, S su i,(P) = a su \ ) (P ) — i H p , understood in a local frame as a 
matrix of first order differential operators. We emphasize the dependence on the 
order of the operator by writing S su b, m (P), so that for P G ’F m (X, £ g) £l^), we 
have •f’sub.m+^-P) = a m (P). 

We shall compute the subprincipal operator of the Laplace-Beltrami operator 
acting on sections of the tensor bundle in Section 4. 

Remark 3.9. For Tb-inner products, the subprincipal operator of an operator P G 
^“(X, £ <g) £7^) acting on £-valued b-half-densities is an element of Diffb( b T*X \ 
0,7 t££), where 7Tb: h T*X \ 0 —>■ X is the projection. In the semiclassical setting, 
P G V%(X,£® fis), we have S suh (P) G Diff^X, ir*£). 

We can now express the symbols of commutators and imaginary parts in a com¬ 
pletely invariant fashion: 

Proposition 3.10. Let P G 4 ,m (X,£®£75) be a TPO with scalar principal symbol 
V- 

(1) Suppose Q G 't m (X, £Tg) £7 2) is an operator acting on £-valued half-densities, 
with principal symbol q. (We do not assume Q is principally scalar.) Then 

a m+ m'-\{P, Q ]) = [S suh (P), q }. 

If Q is elliptic with parametrix Q~, then 

S S ub{QPQ ~) = qSsMP)!- 1 - (3-6) 

(2) Suppose in addition that p is real. Let B be a ^-inner product on £ with 
principal symbol b, then 

P) = Im b S suh (P ), (3.7) 

where lm b S suh (P) = ±(S suh (P) — S su b(P)* b ); we take the adjoint of the 
differential operator S su b{P) with respect to the inner product b on ir*£ and 
the symplectic volume density on T*X. 

Proof. We verify this in a local frame e(x) = {ei(x),..., cn(x)} of £. We compute 

S S ub(P)(^2qjk{x,Quk{x,€)ej(xy\ 

jk 

= E(E Tsnb( P)jk<!k(. i Hp(/:) ai’Cj iqj(H p (ue)ej iqpupCj Hp, 

jt k 

qS SVL b(P) (E ue(x,€)ee{xj S ) 

v t 


while 
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— ^ ^ ^Zjfc^'sub(P)fcl^) Up.&j iqjf\^p('af)^'j iqjpUgCj Hp, 
jl k 

hence S su b(P)q — qS su b(P) = [<7 su b(-P), q] — i^p(q) as an endomorphism (a zeroth 
order differential operator acting on sections of £) of £ in the frame e, which equals 
o-m+m according to the usual (full) symbolic calculus. 

Furthermore, 

Ssub,m(QPQ-) = Ssub, m (P ) + Smb ,m(Q[P,Q~]) 

= Ssub, m (P) + qcJm+m'-l(\P , Q ~]) = S suh , m (P) + q[S suht m(P ), q~ l ] 

= q£sub,m(P)Q , 

noting that Q[P, Q ~] is of order m — 1. 

For the second part, we have S su \ > (P)* b = cr sub (P)* b — (iH p )* b = b~ 1 a su \ 1 (P)*b + 
ib~ 1 (y\ p )*b , where (H p )* is the adjoint of H p as an operator acting on Cf° (T* X \ 0), 
and we equip T*X with the natural symplectic volume density \dxd£\. We have 
(H p )* = — Hp = — Hp since p is real. Therefore, 

S suh (P ) - S suh (P)* b = c S ub(-P) - a sub (P)* b - iHp + ib-'Hpb 
= o- S ub (P) - <r suh (P)* b + ib~ 1 H p {b), 

which indeed gives (3.4) upon division by 2 i. □ 

In particular, (3.7) provides a very elegant point of view for understanding the 
imaginary part of a principally scalar and real (pseudo)differential operator with 
respect to a T-inner product B, as already indicated in the introduction: For 
instance, the principal symbol of the imaginary part Im B P vanishes (or is small 
relative to b = a°(B)) in a subset of phase space if and only if the imaginary part 
of the first order differential operator S su q(P) on T*X \ 0 has vanishing (or small 
with respect to the fiber inner product b of n*£) coefficients in this subset. 

3.4. Interpretation of pseudodifferential inner products in traditional 
terms. We now show how to interpret the imaginary part Im 11 P of an opera¬ 
tor P with respect to a ’F-inner product B in terms of the imaginary part of a 
conjugated version of P with respect to a standard inner product: 

Proposition 3.11. Let B be a 'F -inner product on £. Then for any positive definite 
Hermitian inner product Bq £ C°°(X, Hom(£ ® £1 ?, £ ® 125)) on £, there exists an 
elliptic operator Q £ \F°(A', End(£®f25)) such that B—Q* BqQ £ \F _00 (A, Hom(£® 
£l^,£* ® fis)). 

In particular, denoting by Q~ £ 4 /0 (A, End(£®125)) a parametrix ofQ, we have 
for any P £ 'F m (X ,£ ® f P) with real and scalar principal symbol: 12 

Q(Im B P)Q~ = Im B ° {QPQ~), (3.8) 

and cr m_1 (Im B P) and a " 1-1 (lm B ° (QPQ~)) (which are self-adjoint with respect to 
<r°(B) and Bq, respectively, hence diagonalizable) have the same eigenvalues. 

Proof. In order to shorten the notation, fix a global trivialization of 125 over X 
and use it to identify £ ® Us with £, likewise for all other half-density bundles 
appearing in the statement. Denote the principal symbol of B by b £ «S'° om (T*A' \ 

12 On a symbolic level, this is the same as equation (3.6). 
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0,7r*Hom(£,£ )). We similarly put bo := B$, which is an inner product on tt*£ 
that only depends on the base point. 

We start with on the symbolic level by constructing an elliptic symbol q\ £ 
\ 0, 7T* End(£)) such that b = qlboqi', recall that q\ £ Shom(T*X \ 
0,7T* End(£ )). For t G [0,1], define the Hermitian inner product b t := (l — t)bo + t.b. 
We will construct a differentiable family q t of symbols such that b t = q%boq t for 
t G [0,1]. Observe that for any such family, we have d t b t = b — b 0 = ( d t qt)*b 0 qt + 
qlbod t qt, which suggests requiring d t q t = ) _1 (6 — b 0 ), which we can write 

as a linear expression in q t by noting that ( q J 1 ) -1 = b^qtb^ 1 . Moreover, q 0 = id is a 
valid choice for q t at t = 0. Thus, we are led to define qt, t £ [0,1], as the solution 


of the ODE 


d t q t 


1 

2 


q t b t 1 (b-b 0 ), 


<7o = id • 


Reversing these arguments, for the solution q t we then have ql b^qt = b t for t = 0, 
and both qf boqt and b t are solutions of the same ODE, namely 


dtb t = - ((& - b 0 )b t x b t + b t b t X {b - b 0 )), b 0 = b 0 , 


hence q^b^qt = b t for all t £ [0,1]. 

Let Q i G 'h°(X, End(5)) be a quantization of q±, then we conclude that B — 
QIBqQi £ \H _1 . We iteratively remove this error to obtain a smoothing error: 
Suppose Qk G 4 f0 (X, End(£)) is such that B — Q* k B 0 Q k G \H _fc for some k > 1. 
We will find D k £ a quantization of d k £ S^ k m (T*X \ 0 ,tt*£), such that 

Qk+i ■= Qk + Dk satisfies B — Q£ +1 -BoQfc+i € This is equivalent to the 

equality of symbols 

r k := a~ k (B - Q* k B 0 Q k ) = a~ k (D* k B 0 Q k + Q* k B 0 D k ) = d* k b 0 q 1 + 

which in view of r k = rk is satisfied for dk = j(( b 0 qi)*)- 1 r k ■ We define Q G 
\I/ 0 (X,End(£)) to be the asymptotic limit of the Q k as k -A oo, i.e. Q ~ Qi + 
222—1 D k , which thus satisfies B — Q*B 0 Q G H/ -00 . This proves the first part of 
the proposition. 

For the second part, denote parametrices of B and Q by B~ and Q~ , respectively. 
Then, modulo operators in H/ -00 , we have 


P* B = (. BPB~y = (I Q*B 0 QPQ-Bo 1 {Q-)*y = Q-{QPQ~)* Bo Q , 


hence 

Q(P - P* b )Q~ = ( QPQ~) - ( QPQ~)* Bo 

modulo \1/ —00 . 


□ 


3.5. A simple example. On R" = R Xl x R", , we consider the operator P = 
D Xl + A £ ’L 1 (R ra , C^), where A = A(x,D) £ \E ,0 (R Il ,C Ar ) is independent of xp 
Trivializing the half-density bundle over R n via \dx\?, we can consider P as an 
operator in i[' 1 (R",C JV ® D2). Its principal symbol is oi(P)(x,£) = £1, where we 
use the standard coordinates on T*R", i.e. writing covectors as £ dx, so the Hamilton 
vector field is H CT1 (p) = d x p moreover, in the trivialization of C N by means of its 
standard basis, cr su b(-P)(x, £) = A(x,£). Thus, the subprincipal operator of P is 

S su b(P)(x,0 = A(x,£) - id xx £ Diff 1 (T*R n \ O.tt'C*), 

with A homogeneous of degree 0 in the fiber variables. Suppose we are interested 
in bounding A(p _ p*) on Z := T?,_ 0 ,R" \ 0 relative to a suitably chosen inner 
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product. Let us assume that A(0,£) is nilpotent for all |£| = 1, and that in fact at 
x = 0 and |£| = 1, we can choose a smooth frame ei(£),..., ejy(£) of the bundle 
n*C N — ► T*R n \ 0 so that A(0,£), written in the basis ei(£),..., ejv(£)> is a single 
Jordan block with zeros on the diagonal and ones directly above. Extend the e.j 
by homogeneity (of degree 0) in the fiber variables, and define them to be constant 
in the aq-direction along Z, i.e. ej(x i,0;£) = 6^(0, 0; f), and extend them in an 
arbitrary manner to a neighborhood of Z. 

Now, on Z we have Ae ; j = ej- 1, writing eo := 0. Introduce a new frame e' := e^ej 
with e > 0 fixed, then Ae' } = ee'-_ x . Define the inner product b on n*C N by 
(b(x, £)(e((aq £)), t(e'-(x, £))) = Sij, that is, {e [,..., e' N } is an orthonormal frame 
for b. Then on Z, we find that Im 6 S su \,(P) (which is of order 0) in the frame 
{e[,... ,e' N } is given by the matrix which is zero apart from entries e/2* directly 
above and —e/2 i directly below the diagonal. Thus, defining the ’L-inner product 
B = b(x, D), we have arranged that ||<7o(Im s P)(x, £)||b < e on Z. Since er 0 (Im B P ) 
is self-adjoint with respect to 6, this is really the statement that its eigenvalues are 
bounded from above and below by e and —e, respectively. 

Using Proposition 3.11, we can rephrase this as follows: If Vj denotes the stan¬ 
dard basis of C N and (Bo(vi), t(vj)) = 5^ the standard inner product on 
(the particular choice of an ordinary inner product being irrelevant, see the state¬ 
ment of Proposition 3.11), define the map q(x,£) £ *Sh om (T*M" \ 0,ir*C N ) by 
q(x,£)e'j(: r,£) = Vj. Let Q = q(x,D) and denote by Q~ a parametrix of <5, then 
we find that QPQ~ £ \E ,1 (K”, C N ) satisfies ||cro(Im B ° QPQ~)\\b 0 < e- 

If A has several Jordan blocks not all of which are nilpotent, one can (under the 
assumption of the existence of a smooth family of Jordan bases) similarly construct 
a T-inner product so that the imaginary part of A relative to it is bounded by the 
maximal imaginary part of the eigenvalues of A (plus e) from above, and by the 
minimal imaginary part (minus e) from below. 

4. Subprincipal operators of tensor Laplacians 

Let {M. g ) be a smooth manifold equipped with a metric tensor g of arbitrary 
signature. Denote by TkM = T*M , k > 1, the bundle of (covariant) tensors 

of rank k on M. The metric g induces a metric (which we also call g) on TkM. 
We study the symbolic properties of A k = — trV 2 £ Diff 2 (M, 7/M), the Laplace- 
Beltrami operator on M acting on the bundle TkM. Denote by G £ C°°(T*M) the 
metric function, i.e. G(xT) = ICIgIcc) ’ where G is the dual metric of g. 

Proposition 4.1. The subprincipal operator of Ak is 

S suh (A k )(xT ) = £ Diff 1 (T*M \ 0, n*TkM), (4.1) 

where ' l ~ kM is the pullback connection, with tt: T*M \ 0 —> M being the projec¬ 
tion. 

Proof. Since both sides of (4.1) are invariantly defined, it suffices to prove the 
equality in an arbitrary local coordinate system. At a fixed point xq £ M, introduce 
normal coordinates so that d k gij = 0 at Xq. Then we schematically have 

(AkU)i 1 __i k = gj Ui 1 ,,,i k: jk = T {pk'U J i\...ik,j h ’ du) 

= -g jk dj k Ui 1 ...i k + d(T ■ u) + T • du 
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= -g 3k dj k u il ...i k + T • du + <9r • u, 

with r denoting Christoffel symbols. This suffices to see that the full symbol of A fc 
in the local coordinate system is given by 

c(A k )(x,£) = g jk {x)^j^ k + (x J -x J 0 )e j (x,£)+e(x), 

where ij[x , £) is a linear map in £ with values in End((7fcM) x ), and e(x) is an 
endomorphism of ( T k M ) x . Therefore, cr su b(A fc )(a;o, £) = 0, since d i g :jk (x o) = 0. 
Thus, 

<S' su b(Afc)(a;o, £) = — *H|j|2 = — 2ig jk f k d x j. (4.2) 

We now compute the right hand side of (4.1). First, writing dx 1 = dx 11 g) • • ■<g>dx' lk 
for multiindices I = {i\,... ,ik), we note that sections of Tr*TkM are of the form 
ui(x, f) dx 1 , while pullbacks (under 7r) of sections of T k M are of the form uj(x) dx 1 . 
By definition, the pullback connection V 77 7kM is given by 

Vo TkM (ui(x) dx 1 ) = Vl kM { Ul {x) dx 1 ), Vo TkM (uj(x) dx 1 ) = 0 

^ xd ^xd 

on pulled back sections and extended to sections of the pullback bundle using the 
Leibniz rule; thus, 

v a J kM ( u i(x,Q dx 1 ) = VJ fcM {ui{-,0 dx J ){x), 

xd xd 

Va Sfc TfeM (wj(a;, £) dx 1 ) = %mj( x, £) dx 1 . 

Thus, in normal coordinates at x$ £ M, we simply have V^ ^ M = d x j and 
Vg TkM = d^ k , therefore 

= 2 g j %d xJ 

at Xq, which verifies (4.1) in view of (4.2). □ 

To simplify the study of the pullback connection on n*L~ k M for general k, we 
observe that there is a canonical bundle isomorphism n*TkM = (ff) k tt*T*M; hence 
the connection V 77 7 ~ kM is simply the product connection on (ff) k n*T*M. Therefore, 
if we understand certain properties of 5' su b(Ai), we can easily deduce them for 
S' su b(Afc) for any k. In our application, we will need to choose a positive definite 
pseudodifferential inner product B k = b k {x,D) on the bundle T k M with respect 
to which A k is arbitrarily close to being symmetric in certain subsets of phase 
space. Concretely, this means that we want the operator S su b(Ak) to be (almost) 
symmetric with respect to the inner product b k on ir*7~kM. The following lemma 
shows that it suffices to accomplish this for k = 1: 

Lemma 4.2. Let U C T*M \ 0 be open, and let f £ C°°(U) be real-valued. Fix 
a Hermitian inner product b (antilinear in the second slot) on n*T*M, and define 
R £ End(7 t*T*M) by requiring that 

J (iVu f T M u, v) b da - (u , iV^ 7 M v) b da = J ( u , Rv) b da 

for all u, v £ Cf°(U, tt*T*M), where da is the natural symplectic volume density 
on T* M. There exists a constant C k > 0, independent of U, f and b, such that 
the following holds: If sup^r ||-R||b < e (using b to measure the operator norm of R 
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acting on each fiber) for some e > 0, then the inner product bk = (£) k b induced by 
b on (£) k 7 r*T*M = 7r*7fcM satisfies 

[ (* v h TkM u,v) bk da - [ {u,iV^ TkM v) bk da = f ( u,R k v) bk da , 

Ju Ju Ju 

u,v G C^°(U,ir*TkM), for R k G End(7r*7fcM) satisfying sup^ ||-Rfc||b fc < ke. 


Proof. We show this for k = 2, the proof for general k being entirely analogous. 
Denote S = iV^ T * M , then S 2 = 7 * M acts by ^(uiGt^) = Sui®U 2 +ui®Su 2 - 

Hence using S(au) = aSu + ib\f(a)u for sections u of 7 t*T*M and functions a on 
[7, we calculate 


/ (5 2 (ui <g> u 2 ),v 1 g> v 2 )b 2 da = / {Su 1 ,v 1 ) b {u 2 ,v 2 )b + {u 1 ,v 1 ) b {Su 2 ,v 2 ) b da 
Ju Ju 

= J (ui,S(vi{u2,v 2 )b)^ b + j (u2,S(y 2 {ui,vi)b)^ b da 

+ / (ui <g> u 2 , (R® id + id (g)i?)(ui <g) v 2 ))b 2 da 

Ju 

= / (ui <g> u 2 , S 2 (vi <g> v 2 ))b 2 da - * / H/((ui,ui)b(M 2 ,^ 2 )b) dcr 

Ju Ju 

+ / (ui®U2,R2{vi®V2))boda 

Ju 

= / (ui ® u 2 , <52(ui 0 u 2 ))& 2 + / ( u i ® w 2 ,^ 2 (vi G> u 2 ))b 2 da 

Ju Ju 

with i? 2 = R <8> id + id <g>R, where we used that f (T H/u da = — f n uH/1 da = 0 for 
u G Cf°(U). From the explicit form of i? 2 , we see that ||i?2||b 2 < 2e indeed. □ 


4.1. Warped product spacetimes. Let X be an (n — 1) -dimensional manifold 
equipped with a smooth Riemannian metric h. = h(x,dx), and let a G C°°(X) be 
a positive function. We consider the manifold tf = R t x I, equipped with the 
Lorentzian metric 

g = a 2 dt 2 — h. (4.3) 

On such a spacetime, we have a natural splitting of 1-forms into their tangential 
and normal part relative to a dt, i.e. 

u = ut + ujsiadt. (4.4) 

In this section, we will compute the form of V£| r T M as a 2 x 2 matrix of differential 
operators with respect to this decomposition. For brevity, we will use the notation 
V M := S7 K T M , similarly V' Y := V 71 ’ T x , and we will moreover use the abstract 
index notation, fixing x° = t , and x' = (a; 1 ,..., x n ~ x ) are coordinates on X (inde¬ 
pendent of t ). We let Greek indices run from 0 to n — 1, Latin indices 

i, j, k,... from 1 to n - 1. Moreover, the canonical dual variables 13 £0 =: <7 and 
= (£1 ,... ,£n-i) on the fibers of T*M are indexed by decorated Greek indices 
ft (running from 0 to n — 1) and Latin indices i,j,... (running from 1 to n — 1). 
If an index appears both with and without tilde in one expression, it is summed 


13 Thus, once we discuss Schwarzschild-de Sitter space in the next section, in the region where 
f* = t (which we can in particular arrange near the trapped set), a in the present notation is 
equal to —a in the notation of Section 2. 
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accordingly, for instance ajbj = YHj=i a Xj- Thus, for a section u of n*T*M, we 
have 

Vf u v = Vf u„, Vf u„ = 5^, 

where we interpret V ‘ff as acting on u for fixed values of the fiber variables, i.e. 
viewing a as a family of sections of T*M depending on the fiber variables. As 
before, we denote by G the metric function on T*M , and we let H denote the 
metric function on T* X, interpreted as a (f, cr)-independent function on T* M. 
Lastly, we denote the Christoffel symbols of ( M,g ) by M Tf lv , and those of (X. h) 

i'.v v i'!> 


Lemma 4.3. The Christoffel symbols of M are given by: 

M r o _ n M r o _ -l M r o _ n 

1 00 — u ’ 1 iO — a 1 ij ~ 

M -pfc ^ i kt ^ N1 pfc r» M p/c X ~\^k 

L 00 = oth a e , t i0 = 0, l tj = l y . 


(4.5) 


Proof. We have goo = a 2 , goi = gio = 0 and gij = —hij, and g is ^-independent, 
thus dog^v = 0. Using M V K)lv = \(d^g KV + d v g Kll - d K g l>AJ ), we then compute 

M p _ rj M p _ ^ „ M P _ 

1 000 — u, I OzO — OtOii, 1 Oij — U, 

M p _ _ _ M p _ r» M p _ X p 

k00 QQk') *- kiO '-'j kij ^ kiji 

which immediately gives (4.5). □ 


Proposition 4.4. For the metric g as in (4.3), the subprincipal operator of O i 
(the tensor wave operator acting on 1-forms on M) in the decomposition (4.4) of 
1-forms is given by 

iS sub (Ui)(t,x' ,a, £') 

2a~ 2 adt + cr 2 X _ 2 — X(( H — 2a~ 2 ada 

—2a~ 2 ai\/x a 2a~ 2 ad t + <r 2 H a -2 — H# 




Proof. We start by computing the form of X^u^ and X-fu^ for tangential and 
normal 1-forms. For tangential forms u = dxT with uq = 0, we have 

X(fu 0 = - m T£ 0 u\ = ~a(da,u) H , X^m = d 0 Ui, 

Vfu o = 0, Xfui = X?Ui, Vfu 0 = 0, X% r u i =dp Ui , 


while for normal forms u = dx M with = 0 and Uq = av, we compute 

X^uo = ad t v, Vo I u l = -a t v, 

Xfuo = dj(av) - a jV = adjv, Xf m = 0, X% r u 0 = ad j: v, X~ut = 0. 


Since G = a 2 a 2 
we obtain 


H , we find Hg = 2a 2 ad t + a 2 H a -2 — H h- Using (da, -)h = 


X M — 

y dt ~ 


( ■* 

*V x a 



Moreover, for any / £ C°°(T*X) (we will take / = a -2 and / = iL), viewed as a 
(t, tr)-independent function on T*M , we have Hy = fjdj — fjdj. Hence on tangential 
forms, 
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while on normal forms as above, 

Vh'mo = afjdjV - afjdjv = aH f v, = 0. 

Thus, 



The claim follows. □ 


4.2. Schwarzschild-de Sitter space. We stay in the setting of the previous sec¬ 
tion, and now the spatial metric h has a decomposition 

h = a~ 2 dr 2 + r 2 dw 2 , 

where dui 2 is the round metric on the unit sphere Y = §"^ 2 , with dual metric 
denoted f2; see (2.1). Thus, writing £, resp. r], for the dual variables of r, resp. 
ui £ S n ~ 2 , we have H = a 2 f 2 + r~ 2 \ti\q. Write 1-forms on X as 

u = ut + unO.- 1 dr. (4-6) 

Abbreviate the derivative of a function / with respect to r by /'. Since da = a' dr 
and V x a = a 2 a'd ri we have, in the decomposition (4.6), 

da = ^ a ° a ,^ , iv x a = (0 aa') . 

We will need the Christoffcl symbols of h. We continue using the notation to the 
previous section, except now x 1 = r and = £, while x 2 ,..., x n are r-independent 
coordinates on S™ -2 , and moreover the lower bound for Greek indices is 1, and 2 
for Latin indices. 


Lemma 4.5. The Christoffel symbols of X are given by: 

.V-pl —1 / V1 r\ JV-pl 2 /j 2\ 

r n ol ol , r a o, to [duj 


X -p/c _ i-v X j~\k _ —1 c/c X j~\k _ V -p/c 

^ 11 b, f H ^ i ij f i 


(4.7) 


Proof. We have h\\ = a 2 , /in = hn = 0 and hij = r 2 (dui 2 )ij , and ( doj 2 )ij is 


■ *r 

^-2 ^ _ U.. — O .,.,,1 h . . — JltA,2\.. !A, ,2\ 

r-independent. We then compute 

- Y T m = —a~ 3 a', X T U1 = 0, ' Y Tijj = -r(du 2 )ij, 

X r fc n = 0, x T kil = r(duj 2 ) k i, X T kij = r 2Y T kij , 

which immediately gives (4.7). □ 

We are only interested in the subprincipal operator of llli at the trapped set, 
which we recall from (2.6) to be the set 

r = {r = r p , £ = 0, a 2 = 4' 2 |r 7 | 2 }, where = or -1 ,4/'(r p ) = 0. (4.8) 

Thus, at T, we have 

H# = 2a 2 fd r - 2aa'^ 2 d^ + 2 r _3 |? 7 | 2 d£ + r ~ 2 \^\ v \2 = 2W 3 |?y| 2 <9 5 + r~ 2 H\ v ^, 

while cr 2 H a -2 = 2a 2 a~ 3 a'd^. Now a~ x a' = (r\E') _1 (r\E') / = r _1 at r = r p , therefore 
a 2 a~ 3 a! = r~ 3 \r/\ 2 , and we thus obtain 

cr 2 H a -2 — H h = — r ~ 2 H| r( | 2 at r. 

Notice that |?y| 2 £ C°°[T*Y) is independent of (r, £). 


(4.9) 
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Lemma 4.6. For a function f E C°°(T*Y), viewed as an (r, £)-independent func¬ 
tion on X, we have 

Vj!C ar(i Hf duj 2 )\ 
ar _1 i H/ H f ) 

in the decomposition (4.6) of 1-forms on X. 

Proof. On tangential forms u, i.e. ui = 0, we have 

V' Y rti = —r~ 1 Uj, Vf Ui = V] Ui , X~u i = 0, V~ Y iii = djUi, 

thus using H/ = f~dj - fjd~ p we get 14 

V mi = -r _1 /jitj = -r _1 u(H/) = -r~ 1 i Hf u, = V^iq. 

On normal forms it, i.e. iti = a _ 1 u, iq = 0, we compute 

V^iti = a~ 1 djV, V‘ Y it* = ra(duj 2 )ijV , = a” 4 cLi;, Vj'iq = 0, 

hence 

V^ui = a -1 fjdjV - a~ 1 f j djV = a _1 H/t;, 

Vtf = fjra(du} 2 )ijV = ar(i Hf duj 2 )v. 

The claim follows immediately. □ 

Combining Proposition 4.4 and Lemma 4.6, we can thus compute the subprinci¬ 
pal operator of Di acting on 1-forms (sections of the pullback of T*M to T*M \ 0) 
decomposed as 

u = utt + utnoT 1 dr + u^oidt. (4-10) 

In view of (4.9), we merely need to apply Lemma 4.6 to f = \r]\ 2 , in which case 
H/ = 2 fF fe ? 7 i dfc - diCl^rjjTjkdj, so iHj = 2 i v on 1 -forms (identifying the 1 -form rj 
with a tangent vector using the metric dw 2 ), while iy\ f dijj 2 = 2 77 . Thus, we obtain: 

Proposition 4.7. In the decomposition (4.10), the subprincipal operator ofO 1 on 
Schwarzschild-de Sitter space at the trapped set T is given by 

iSsubif-^l) 

( 2 a~ 2 odt — r~ 2 Vu — 2 ar~ 1 n 0 

h m 2 ' 

— 2 ar 3 i v 2a 2 odt~r 2 Hi^p —2r l o 

y 0 —2r _1 cr 2a~ 2 odt — r _2 H |, ; |2 

(4.11) 

Since Di is symmetric with respect to the natural inner product G on the 1- 
forrn bundle, which in the decomposition (4.10) is an orthogonal direct sum of 
inner products, G = (—r _ 2 fl) ® (—1) ® 1, the operator 5 su b(ni) is a symmetric 
operator acting on sections of n*T*M over T*M \ 0 if we equip ir*T*M with the 
fiber inner product G and use the symplectic volume density on T*M \ 0. 

The matrix — 2r~ 2 s, with 

/ 0 ihr 2 ?7 0 \ 

s = I —0 ro I , 
y 0 ro 0 J 

44 We use that 7 t*T*X can be canonically identified with the horizontal subbundle of T* (T*X). 
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of 0-th order terms of 5 su b(ni) is nilpotent, which suggests in analogy to the 
discussion in Section 3.5 that the imaginary part of SsubG^i) with respect to a 
Riemannian fiber inner product can be made arbitrarily small. Indeed, for any 
fixed e > 0, define the ‘change of basis matrix’ 

( id 

q= 0 

\—e~ 2 \i]\~ 1 'l’ 2 r 2 i ri 


0 0 
e -!^ r 2 o 

0 e _2 |ry| _1 4'r 3 cr 


then 

/0 eq 0 \ 

qsq- 1 = 0 0 e\rj\ . 

\0 0 0 / 

In order to compute gS'sut^Di)? -1 , we note that the diagonal matrix of f-derivatives 
in (4.11) commutes with q , and it remains to study the derivatives along more 
specifically, q has a block structure, with the columns and rows 1, 3 being the first 
block and the (2, 2) entry the second, and the (2, 2) block is an ^-independent 
multiple of the identity, hence commutes with the relevant (2, 2) entry ir~ 2 of 
SUbCHi). For the 1,3 block, we compute 


V h„ 


H 


0 

M 2 


id 

’ l —e” 2 |?7|” 1 \I ,2 r 2 ^ 


2 |??| 1 4'r 3 CT 




0 0 N 
H|.^|2i^ 0 


(4.12) 


= e 2 ^ 2 r 2 \rj\ 1 

Now Vh and H^p are the restrictions of the pullback connection Vof 
the full form bundle to 1-forms and functions, respectively, and the latter commutes 
with i r] . since by Proposition 3.10, 


0 = s sub ([□,£]) = -i[s sub (□),*„] = -[vs^f 

where □ denotes the Hodge d’Alembertian on the form bundle and 5 is the codif¬ 
ferential. Thus, (4.12) in fact vanishes, and therefore 

qS^n^q- 1 

2 a~ 2 ad t — r _2 V?, —2 r 2 eq 0 

i”i „ ' n 

0 2a 2 ad t —r 2 H|.^ 12 — 2 r 2 e|? 7 | 

0 0 2a~ 2 ad t — r _2 H| J) p 




Equip the 1-form bundle over Ad in the decomposition (4.10) with the Hermitian 
inner product 

B 0 = n®l®l, (4.13) 

then qSsubiD^q^ 1 has imaginary part (with respect to Bq) of size O(e). Put 
differently, 5 , su b(Di) has imaginary part of size 0(e) relative to the Hermitian inner 
product b := Bo(q-,q-), which is the symbol of a pseudodifferential inner product 
on n*T*Ad. We can now invoke Lemma 4.2 on a neighborhood of T D {|er| = 1} and 
use the homogeneity of q , b and *S' su b(Di) to obtain: 


Theorem 4.8. For any e > 0, there exists a (positive definite) i *-independent 
pseudodifferential inner product B = b(x,D) on 7 (thus, b is an inner product 
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on n*T k M, homogeneous of degree 0 with respect to dilations in the base T*M\0), 
such that 


sudIctI” 1 -(S^(a,,)~ S^(Di.)* b ) < e, 



where T is the trapped set (4.8). Put differently, there is an elliptic 'ItDO Q, in¬ 
variant under f* -translations, acting on sections ofT k M, with parametrix Q~, such 
that relative to the ordinary positive definite inner product (4.13), we have 



< e. 


By restriction, the analogous statements are true for □ acting on subbundles of the 
tensor bundle on M , for instance differential forms of all degrees and symmetric 
2-tensors. 

By the ^-translation invariance of the involved symbols, inner products and 
operators, this is really a statement about U/b-inner products, and Q is a b-pseudo- 
differential operator; see the discussion preceding Theorem 2.1 for the relationship 
of the stationary and the b-picture. 

Remark 4.9. Adding a 0-th order term to □ does not change □ or its imaginary 
part at the principal symbol level, thus does not affect the subprincipal operator of 
□ either; therefore, Theorem 4.8 holds in this case as well. 

Adding a first order operator L (acting on sections of T k M), which we assume to 
be f-independent for simplicity, does affect the subprincipal operator, more specif¬ 
ically its 0-th order part, since S' su b(n + L) = S S ub(D) + or (A). Thus, if or (A) 
is small at T, we can use the same dt-inner product as for □ and obtain a bound 
on Im b «S' su b(n + A) which is small, but no longer arbitrarily small. However, the 
bound merely needs to be smaller than i n /2, see (2.5), which does hold for small 
A. 

If we do not restrict the size of A, we can still obtain a spectral gap, provided 
one can choose a Ik-inner product as in Theorem 4.8, again with e > 0 sufficiently 
(but not necessarily arbitrarily) small. This is the case if the 0-th order part of 
S mb (D + A) is nilpotent (or has small eigenvalues) and can be conjugated in a t- 
independent manner to an operator which is sufficiently close to being symmetric, 
in the sense that it satisfies the bound (2.5) with □ replaced by □ + A. 

We remark that the subprincipal operator iiS su b(D) = + *cr su b(G) induces 

a notion of parallel transport on n*T k M along the Hamilton flow of H^. As a 
consequence of the nilpotent structure of S' su b(n) at the trapped set, parallel sec¬ 
tions along the trapped set grow only polynomially in size (with respect to a fixed 
A invariant positive definite inner product), rather than exponentially. Parallel sec¬ 
tions as induced by 5' su b(D + A), with A as in Remark 4.9, may grow exponentially, 
with their size bounded by C'e K l cr l* for some constants C > 0 and k, where the 
additional factor of \a\ in the exponent accounts for the homogeneity of the parallel 
transport. If such a bound does not hold for any n < v m i n /2, the dispersion of 
waves concentrated at the trapped set caused by the normally hyperbolic nature 
of the trapping is expected to be too weak to counteract the exponential growth 
caused by the subprincipal part of □ + A, and correspondingly one does not expect 
a spectral gap. Notice that the growth of parallel sections is an averaged condition 
in that it involves the behavior of the parallel transport for large times, while the 
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choice of 'JZ-inner products as explained above is a local condition and depends 
on the pointwise structure of S' su b(D); thus, establishing spectral gaps only using 
averaged data is an interesting open problem, even in the scalar setting. 
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